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ABSTRACT 

Short lived radionuclides (SLRs) like 26 Al are synthesized by massive stars and are a byproduct of star 
formation. The abundances of SLRs in the gas of a star-forming galaxy is inversely proportional to its gas 
consumption time. The rapid evolution of specific star formation rate (SSFR) of normal galaxies implies they 
had mean SLR abundances ~ 3 — 10 times higher at z = 2. During the epoch of Solar System formation, the 
background SLR abundances of the Galaxy were up to twice as high as at present, if SLR yields from massive 
stars do not depend on metallicity. If SLRs are homogenized in the gas of galaxies, the high SSFRs of normal 
galaxies can partly explain the elevated abundance of SLRs like 60 Fe and 26 Al in the early Solar System. 
Starburst galaxies have much higher SSFRs still, and have enormous mean abundances of ~ D A1 ( zo Al/-'Al 
« 10~ 3 for Solar metallicity gas). The main uncertainty is whether the SLRs are mixed with the star-forming 
molecular gas: they could be trapped in hot gas and decay before entering the colder phases, or be blown 
out by starburst winds. I consider how variability in star-formation rate affects the SLR abundances, and 
I discuss how SLR transport may differ in these galaxies. The enhanced 26 Al of starbursts might maintain 
moderate ionization rates (10~ 18 - 10~ 17 sec" 1 ), possibly dominating ionization in dense clouds not penetrated 
by cosmic rays. Similar ionization rates would be maintained in protoplanetary disks of starbursts, if the SLRs 
are well-mixed, and the radiogenic heating of planetesimals would likewise be much higher. In this way, galaxy 
evolution can affect the geological history of planetary systems. 



Subject headings: galaxies:starburst — 
thesis, abundances ■ 



1. INTRODUCTION 



galaxiesTSM — galaxies: evolution — nuclear reactions, nucleosyn- 
- solar system: formation — stars: formation 



Massive, young stars are sites of nucleosynthesis, not just 
of the stable nuclides, but radioactive isotopes as well. Long- 
lived radioisotopes with decay times of billions of years or 
longer, such as 40 K, mix into the interstellar medium (ISM), 
accumulating and pr ovide a very low level source of energetic 
radiat ion in all gas (lCameronlll962t lUmebavashi & Nakanol 
119811) . More unstable isotopes are also synthesized; some 
have decay times of a few years or less and cannot reach most 
of the ISM. In between those two extremes are the short lived 
radionuclides (SLRs) of the ISM: with <~ Myr decay times, 
they can be present in a galaxy's gas but only as long as 
star formation replenishes their abundances. The most promi- 
nent of the SLRs is 26 Al, which is detected in gamma-ray de- 
cay li nes from star-format i on regions through out the Milky 
Way dMahonev et al.|[l98i iDiehl et alj|2006al) . Another SLR 
that has recently b een detected in gamma-ray lines is 60 Fe 
dHarris et al.ll2005b . 

SLRs are not just passive inhabitants of the ISM. By releas- 
ing energetic particles and radiation when they decay, they 
inject power that can heat or ionize the surrounding matter. In 
the Milky Way's molecular clouds, the radioactivity is over- 
whelmed by that of cosmic rays, which sustain an ionization 
rate of £cr ~ 5 x 10~ 17 sec" 1 . Rapidly star-forming galaxies 
known as starbursts may have elevated levels of cosmic rays 
and ionizati on rates a thousand times or more higher t han the 
Milky Way ( ISuchkov et alJ[l993l:lPapadopoulosl2010l) . How- 
ever, it is possible that cosmic rays are stopped by the abun- 
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dant gas in starbursts before they enter the densest molecu- 
lar gas. Gamma rays can provide ionization through large 
columns; but while the gamma-ray ionization rate can reach 
up to C7 ~ 10~ 16 sec" 1 in the densest starbursts, in mo st star- 
bursts gamma rays sustain relatively weak ionization (lLackil 
120121) . SLRs like 26 Al can in principle provide ionization 
through any column of gas, and if abundant enough, main- 
tain moderate levels of ionization. The major open question 
for SLR-induced ionization is how well mixed the SLRs are 
with the gas of the galaxy, a process which takes anything 
up to 100 Myr; if the mixing times are longer than a few 
Myr, the SLRs are not abundant in most star-forming cores 
dMever & Clavtonll2000l iHuss et al.ll2009l) . 

Meteorites recording the composition of the primordial 
Solar System demonstrate that SLRs were present dur- 
ing its formation. Assuming the SLRs were not cre- 
ated in situ by energetic radiation from the Sun dLee et al.l 
1998), the SLRs provide evidence that the Solar Sys- 
tem formed near a star-forming region with young mas- 
sive stars (e.g.. lAdamsll2010l) . In fact, 26 Al was overabun- 
dant by a factor ~ 6 in the primordial Solar System, with 
X( 26 A1) sa 10" 10 ( 26 A1/ 27 A1 w 5 x 10" 5 ), compa red to its 
present day abundances in th e Milky Way (e.g . . iLee et al 
19771 iMacPherson et al.|[l99l IDiehl eTal]l2006ai; fHuss et al 



2009). Their quick decay time also indicate the Solar Sys- 
tem formed quickly, within about a few Myrs. SLRs, 
particularly 26 Al, were a pr i mary source of ionization in 
the Solar Nebula dStepinskil [19921: iFinocchi & Ga5l [l997l 
lUmebavashi & Nakanol 20091) . affecting the conductivity and 
ultimately accretion rate in the protoplanetary disk. More- 
over, 26 Al and other SLRs may have regulated the early geo- 
logical evolution of the Solar System by being a major source 
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of heat in early planetesimals , driving their differentiation 
and rock metamorphism (e.g.. iHutcheon & Hutchisonlll989t 
iGrimm & McSween|[T99llShukolvukov & Lugmairill993l) . 

The contemporary Milky Way, with a star-formation rate 
(SFR) of a few stars per year is not a typical environment 
for most of the star-formation in the history of the Universe, 
however. Roughly 5-20% of star-formation at all times oc- 
curred in rapid starburst s mainly driven by galaxy-galaxy in- 
teract ions and mergers dRodighiero et al.ll201 ll ISargent et alj 
120121) . Furthermore, most of the star-formation in "normal" 
galaxies occurred in massive galaxies with a much higher star- 
formation rate (> 10 M Q yr" 1 ) at redshifts z of 1 and higher, 
when most star-formation took place (e.g., iMagnelli et al.l 
120091) . These high star-formation rates translate into large 
masses of SLRs present in these galaxies. I shall show that 
26 Al in these galaxies, if it is well mixed with the gas, can 
sustain rather high ionization rates in their ISMs. This has 
consequences for both star formation and planet formation. 

When necessary, I assume a Hubble constant of Ho = 
72 km s" 1 Mpc -1 , a matter density of VLm = 0.25, and a cos- 
mological constant SI a = 0.75 for the cosmology. 

2. THE EQUILIBRIUM ABUNDANCE OF SLRs 

In a one-zone model of a galaxy, which disregards spa- 
tial inhomogeneities, the complete equation for the SLR mass 
M SL r in the ISM is 

dMsLR n _ M SLR (r) 



dt 



Gslr(0- 



TSLR 



(1) 



where tslr is the lifetime of the SLR in the galaxy. Qslr(?X 
the injection rate of the SLR, depends on the past star- 
formation history: 



GsLr(0 = 



7 S LR(f-f')xSFR(fVf'- 



(2) 



For a coeval stellar population of age t, the yield Fslr(0 is the 
mass ejection rat e of the SLR into the interstellar medium per 
unit stellar mass dCervifio et alj2000b . 

If there are no big fluctuations in the star-formation rate 
over the past few Myr, then the SLR abundance approaches 
a steady-state. The equilibrium mass of a SLR in a galaxy is 
proportional to its star-formation (or supernova) rate averaged 
over the previous few Myr. We can then parameterize the 
injection of SLRs in the ISM by a yield Tslr per supernova: 



T SL r = e I Y(t")dt", 
lo 



(3) 



regardless of whether SNe are actually the source of SLRs. 
The e factor is the ratio of the supernova rate Tsn and star- 
form ation rate. Then th e equilibrium SLR mass is given by 
(e.g., iDiehl et alj|2006al) 



TsnTslrtslr- 



(4) 



The supernova rate is proportional to the star-formation rate 
(SFR), so Tsn = eSFR. The abundance of an SLR is given 
by X SL r = Ms LR m // /(M gas m S LR), where m SL R is the mass of 
one atom of the SLR and M sas is the gas mass in the galaxy. 
Therefore the abundance of an SLR is 



^SLR = £ 



SFR T SLR T SLR m H 



M, 



gas 



?«SLR 



(5) 



The quantity M gas /SFR = r gas is the gas consumption time. 
Note that it is related to the specific star formation rate, 



SSFR = M./SFR, as r gas = / gas /((l -/ gas )SSFR), where M* 
is the stellar mass and / gas =M gas /(M gas +M + ) is the gas frac- 
tion. Therefore, we can express the equilibrium mass of the 
SLR in a galaxy as 

£ T SLR ''"SLR m H l _ /gas„„™, TsLRTSLR^H 

Xslr = = £ 7 SSFR (6) 



Tgas'«SLR 



'«SLR 



Finally, the ratio of SLR abundance in a galaxy to that in the 
Milky Way is 

v T MW 

a slr T gas tslr 



yMW 
SLR 



gas 
'gas SLR 

\-f f 

1 ./gas J i 



MW 



SSFR 



/gas SSFR MW r MW ' 



(7) 



with a MW superscript referring to values in the present day 
Milky Way. Thus, galaxies with short gas consumption times 
(and generally those with high SSFRs) should have high abun- 
dances of SLRs. The reason is that in such galaxies, more of 
the gas is converted into stars and SLRs within the residence 
time of an SLR. 

The greatest uncertainty in these abundances is the resi- 
dence time tslr- In the Milky Way, these times are just the 
radioactive decay times, defined here as the e-folding time. In 
starburst galaxies, however, much of the volume is occupied 
by a hot, low density gas which forms into a galactic wind 
with characte ristic speeds v of several hundred kilometers per 
second (e.g.. iChevalier & Clegd [l985T: iHeckman etai1[T990T: 
iStrickland & Heckmanll2009l) . If massive stars emit SLRs at 
random locations in the starburst, most of them will dump 
their SLRs into the wind phase of the ISM. The wind-crossing 
time is T w j n d = 330kyr(/!/100pc)(v/300km s" 1 )" 1 , where h is 
the gas scale-height. The equilibrium time in starburst galax- 
ies is then r = [t^^+t'I^]' 1 . Furthermore, the SLRs ejected 
into the wind may never mix with the molecular gas, so the 
fraction of SLRs injected into the molecular medium may be 
C 1 (I discuss this issue further in section 13. 2\ . However, 
very massive stars are found close to their birth environments 
where there is a lot of molecular gas to enrich, and these may 
be the source of 26 Al, as supported by the correlation of the 
1.809 MeV 26 Al decay li ne emission and fr ee-free emission 
from massive young stars (Knodlseder 1999b . 

Turning to the specific example of 26 Al, I note that the 
yield of 26 Al i s thought to be T ai- 26 « 1.4 X 10" 4 M 
per supernova (IDiehl et al.l l2006al) . For a Salpeter ini- 
tial mass function from 0.1 - 100 M Q , the supernova 
rate is r SN = 0.0064 yr _1 (SFR/M Q yr" 1 ), or r SN = 
0.11 yf-'CLnR/lO 11 L ) in terms of the t otal infrared (8 - 
1000 /xtn) luminosity Ltir of starbur sts dKennicuttl 119981 : 
iThompson. Ouataert. & Waxmari 120071) . If I suppose all of 
the 26 Al is retained by the molecular gas, so that the residence 
time is the 26 Al decay time of 1 .04 Myr, then the equilibrium 
abundance of 26 Al in a galaxy is just 



X( 2b Al) = 3.4x 10" 



^L) =1.7x10" 
Gyry 



20 Myr 



(8) 



2.1. High-Redshift Normal Galaxies 

The star-formation rates and stellar masses of normal star- 
forming galaxies lie on a "main sequence" with a charac- 
teristic SSFR that varies weakly, if at all, with stellar mass 
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(e.g.. iBrinchmann et al.l [2004b . However, the characteristic 
SSFR evolves rapidly with redshift ( e.g., iDaddi et al.l [20071: 
iNoeske et all2007t iKarim et alJl20TTh. with SSFR oc (1 + z) 2 8 
out to z « 2.5 - a rise of factor ~ 30 dSargent et al.ll2012l) . At 
z > 2.5, t he SSFR of the main sequence then seems to remain 
constant ([Gonzalez et al.ll20lol) . 

Countering this rise in the SSFR, the gas fractions of nor- 
mal galaxies at high z were also higher: the high equilib- 
rium m asses of SLRs a r e dilu ted to some extent by higher gas 
masses. IHopkins et al.l (120101) p rovide a conveni ent equation, 
motivated by the Schmidt law dKennicuttll 19981) . to describe 
the evolution of gas fraction: 

/gas(z) = / [l-fc(z)Ao)(l -/o /2 )]" 2/3 , (9) 

assuming a gas fraction /o at z = 0, with a look back time 

of t L (z) = ft dz' I [H ( 1 + z') \/n a + 1 + z') 3 l a nd a current 
cosmic age of to (see also IHopkins et all 120091) . Since the 
gas fractions of normal galaxies at present are small, the 
evolution at low redshifts can be approximated as / gas (z) = 

/ [i-te(z)Ao)r 2/3 . 

After calculating the mean abundances of SLRs in normal 
galaxies, I find that the rapid SSFR evolution overwhelms the 
modest evolution in / gas at high z: the SLR abundances of nor- 
mal galaxies evolves quickly. These enhancements are plotted 
in Figure Q] 

Observational studies of high redshift main sequence galax- 
ies indicate a slower evolution in r gas , resulting from a quicker 
evolution of / gas . Although equation [9] implies that / gas was 
about twice as high ten billion years ago at z ~ 2, massive disk 
galaxies are observed with gas fractions of ~ 40 -50%, which 
is 3 to 10 times greater than at presen t (e.g.. iTacconi et al.l 
I2010t IDaddi et al.ll2010l) . According to lGenzel et al.l (120101) . 
the typical (molecular) gas consumption time at redshifts 1 to 
2.5 was ~ 500 Myr. In the IDaddi et "all (1201 Oh sample of BzK 
galaxies at z ~ 1 .5, gas consumption times are likewise ~ 300 
- 700 Myr. To compare, the molecular gas consu mption times 
at the present are estimated to be 1 .5 to 3 Gyr (iDiehl et al. 
2006at iGenzel et all 120 lot iBigiel et alj|201 U iRahman et al" 
20121) . implying an enhancement of a factor 3 to 6 in SLR 
abundances at z > 1 . But note that the BzK galaxies are not 
the direct ancestors of galaxies like the present Milky Way, 
which are less massive. The SSFR, when observed to have 
an y mass dependence , is greater in low mass galaxies at all 
z jSargent et al.ll2012b . This means that lower mass galax- 
ies at all z have shorter T pas , as indica ted by observations of 
present galaxies (ISaintonge et al.l201 ll) . The early Milky Way 
therefore may have had a gas consumption time smaller than 
500 Myr. 

So far, I have ignored possible metallicity Z dependencies 
in the yield T of SLRs. It may be generally expected that 
star-forming galaxies had lower metallicity in the past, since 
less of the gas has been processed by stellar nucleosynthe- 
sis. However, observations of the age-metallicity relation of 
G dwarfs near the Sun reveal that they have n early the same 
metallicity at ages approaching 10 Gy r (e.g., lTwarog|[l980l 
lHavwoodfeooa iHolmberg et al.ll2007l) . though the real sig- 
nificance of the lack of a trend remains unclear, since there 
is a wide scatter in metallicity with age (see the discussion 
by iPrantzosi 120091) . Observations of external star-forming 
galaxies find weak evolution at constant stellar mass, with 
metallicity Z decreasing by ~ 0.1 -0.2 dex pe r unit redshift 
dLillv et al.ll200llKobulnickv & Kewle\ll2004T) . 

After adopting a metallicity dependence of Z(z) = Z(0) x 




z 



FIG . 1 . — Plot of the SLR abundance enhancements in normal galaxies lying 
on the "main sequence", for a gas fraction evolution described by equation|9j 
The rapid evolution of SSFRs leads to big enhancements of SLRs at high z- 
Even during the epoch of Solar System formation, the mean SLR abundance 
was twice the present value. The different lines are for different /gas at z = 
0, assuming SLR yields are independent of metallicity: 0.05 (dotted), 0.1 
(solid), 0.2 (dashed). The shading shows the abundances for 0.05 < / gas < 
0.2 when the SLR yield depends on metallicity, assuming a 0.2 dex decrease 
in metallicity per unit redshift. 

\0~°- 2z (0.2 dex decrease per unit redshift), I show in the re- 
vised SLR abundances Figure Q] assuming that the SLR yield 
goes as Z" 1 , Z" 5 , Z 5 , Z, Z , and Z 2 . If the yields are 
smaller at lower metallicity, the SLR abundances are still ele- 
vated at high redshift, though by not as much for metallicity - 
independent yields. As an example, the yield of 26 Al in 
the winds of Wolf-Rayet stars is believed to scale as T oc 
Z 15 (lPalaciosetalJl2005l) . According to iLimongi & Chieffil 
d2006l) . these stellar winds contribute only a minority of the 
26 Al yield, so it is unclear how the 26 Al yield really scales. 
iMartin et ail ( 120101) considered the 26 Al and 60 Fe yields from 
stars with half Solar metallicity. They found that, because 
reduced metallicity lowers wind losses, more SLRs are pro- 
duced in supernovae. This mostly compensates for the re- 
duced wind 26 Al yield, and actually raises the synthesized 
amount of 60 Fe. 

2.2. Starbursts 

The true starbursts, driven by galaxy mergers and galaxies 
interacting with each ot her, represent about ~ 10% of star for- 
matio n at all redshifts dRodighiero et al. 2011; Sarge ntet alj 
120121) . They have SSFRs that are up to an order of mag- 
nitude higher than z = 2 normal galaxies. The mean, back- 
ground abundances of SLRs in starbursts are therefore about 
100 times greater than the present day Milky Way. 

I show the 26 Al abundances in some nearby starburst galax- 
ies in Table Q] In the Galactic Center region, the 26 Al abun- 
dance is only twice that of the present Milky Way as a whole. 
However, the 26 Al abundances are extremely high in the other 
starbursts, ~2x 10~ 9 , about twenty times that of the primor- 
dial Solar nebula. The 26 A1/ 27 A1 ratio in these starbursts is 
also very high. Assuming Solar metallicity with an 27 Al abun- 
dance of log 10 [jV( 27 Al)/M7/)] = -5.6 (IDiehl et alJ l2006al) . 
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TABLE 1 

6 Al Abundances and Associated Ionization Rates 



Starburst 


SFR 




M ct| 
m M-26 


M h 


T gas X( 26 Air 


( 26 A1 / 27 Al) b 


fe c 
^Al-26 


f*~t d 




(M Q yr" 1 ) 


(yr 1 ) 


(M Q ) 


(M s ) 


(Myr) 




(sec - ) 


(sec -1 ) 



Milky Way (z = 0) c 


3.0 


0.019 


2.8 


4.5 x 10' 


1500 


2.4 X 10-" 


9.4 x 10- 6 


1.9 x 10- 20 


7.1 x 10- 20 


Galactic Center CMZ f 


0.071 


4.6 x 10- 4 


0.067 


3 x 10 7 


420 


8.6 x 10-" 


3.4 x 10- 5 


6.9 x 10- 20 


2.6 x 10- 19 


NGC 253 Core°- h 


3.6 


0.023 


3.3 


3 x 10 7 


8.3 


4.3 x 10- 9 


1.8 x 10- 3 


3.4 x 10-' 8 


1.3 x 10 -17 


M828- 1 


10.5 


0.067 


9.8 


2x 10 s 


19 


1.9 x 10-' 


7.5 x 10- 4 


1.5 x 10- 18 


5.7 x ur 18 


Arp 220 Nuclei 


50 


0.3 


44 


10" 


20 


1.7 x 10- 9 


6.7 x 10- 4 


1.3 x 10-' 8 


5.0 x 10 -18 


Submillimeter galaxy k 


1000 


6.4 


930 


2.5 x 10 10 


25 


1.4 x 10- 9 


5.7 x 10- 4 


1.1 x 10- 18 


4.3 x 10 -18 


BzK Galaxies 1 


200 


1 


200 


7x 10 10 


400 


1 x 10 -10 


4 x 10- 5 


8 x 10- 20 


3 x 10- 19 



a Mean abundance of 26 Al, calculated assuming the 26 Al is well-mixed with the gas and resides there for a full decay time (instead of, for example, 
a wind-crossing time). 

Calculated assuming Solar metallicity with log 10 [W( 27 Al)/A'(//)] = -5.6. 
c Ionization rate from 26 Al with the derived abundance, with ionization only from MeV positrons released by the decay, assuming effective 
stopping. 

d Ionization rate from 26 Al with the derived abundance, where ionization from the 1.809 MeV decay line and 0.511 keV positron annihilation 
gamma rays is included, assuming t hey are all s t opped. 

"■ Supernova rate and gas mass from Diehl et al. (2006a); SFR calculated from supernova rate using Salpeter IMF for consistency. 



f Inner 100 pc of the Milky Way. SFR and Tsn from IR luminosity in Launhardt et al. (2002); gas mass from Molinari et al. (2011) 
Pierce-Price et al. (2000) gives a gas mass of 5 X 10 7 Mq. 
g SFR and Tsn from IR luminosity in Sanders et all J2003I) . 
h Gas mass fromlHarrison et al. ( 19991). 

' Gas mass from lWeifi et al.H2Q01l) „^__^ 

j Assumes IR luminosity of 3 X 10" Lq for SFR and Tsn and gas mass given in Downes & Solomon 1 1998). 

Typical gas mass and SFR of submillimeter galaxies from Tacconi et al. (20(H). 
' Mean SFR and gas mass of the 6 BzK galaxies in Daddi et al. 1 201C), which are representative of main sequence galaxies at z ~ 1-5. 



this ratio is ~ (0.6- 1.8) x 10 3 . Again, this ratio for Solar 
metallicity gas is ~ 10-30 times higher than that of the early 
Solar Nebula, - 5 x 10" 5 . 

3. SYSTEMATIC UNCERTAINTIES 

3.1. Effects of Variable Star-Formation Rates 

The steady-state abundance (eqn. |4]i is only appropriate 
when the star-formation rate is slowly varying on timescales 
of a few Myr. Since young stellar populations produce SLRs 
for several Myr, and since 26 Al and 60 Fe themselves survive 
for > 1 Myr, the injection rate of SLRs is smoothed over those 
timescales (equation [2]). Very high frequency fluctuations in 
the SFR therefore have little effect on the abundance of SLRs. 
In the opposite extreme, when the fluctuations in SLRs are 
slow compared to a few Myr, we can simply take the present 
SFR and use it in equation [4] for accurate results. However, 
intermediate frequency variability invalidates the use of equa- 
tion|4] and can result in the SLR abundance being out of phase 
with the SFR. 

Normal main sequence galaxies at high redshift built up 
their stellar population s over Gyr times, evolving secularly 
(c.f.. lWuvts etaDIIoTl . They are also large enough to con- 
tain many stellar clusters, so that stochastic effects average 
out. It is reasonable to suppose that they have roughly con- 
stant SFRs over the relevant timescales. True starbursts, on 
the other hand, are violent events that last no more than 
^100 Myr, as evinced by their short T gas . They are rela- 
tively small, so stocha stic fluctuations in their star-fo rmation 
rates are more likely. IForster Schreiber et al.l (l2003t studied 
the nearby, bright starburst M82 and concluded that its star- 
formation history is in fact bursty. The star-f ormation his- 
tories of other starbursts are poorly known, but iMavva et al.l 
( 120041) present evidence for large fluctuations on ~ 4 Myr 
times. 

I estimate the magnitude of these fluctuations for the proto- 
typical starburst M82 with the full equation for SLR mass in 
a one-zone model (equationQ]). The solution to equationQ]for 



Mslr is 



where 



M SL r(/)= f SFR(t')xm SLR (t')dt', (10) 

J -oo 



m SLR (t')= / T S LR(f")exp 



t'-t" 

7SLR 



dt". (11) 



The quantity msLR(f') represents the SLR mass in the ISM 
from a coeval stellar population of unit mass and age t 1 . It is 
given bv lCervino etail (120001) and lVoss et all (120091) for 26 Al 
and 60 Fe. 

I use the star -f ormat ion history derived by 
IForster Schreiber etal.l (12003b for the "3D region" of 
M82, which consists of two peaks at 4.7 Myr ago and 8.9 
Myr ago. The p eaks are modelled as Gaussian s with the 
widths given in IForster Schreiber et al.l (120031) (standard 
deviations a of 0.561 Myr for the more recent burst, and 
0.867 Myr for the earlier burst). I convert the star-formation 
rate from a Salpeter I MF f rom 1 to 100 M Q given in 
IForster Schreiber et al.l (120031) to a Salpeter IMF from 0.1 
to 100 M Q for consistency with the rest of the paperQ This 
region does not include the entire starburst; it has roughly 1 /3 
of the luminosity of the starburst, but the stellar mass formed 
within the 3D region over the past 10 Myr gives an average 
SFR of 10 Mq y r -1 in the IForster Schreiber et al.l (12003b 
history. Note that iRodriguez-Merino et al 1 1201 3) derives a 
different age distrib u tions for stellar clusters (comp are with 
ISatvapal et aDl997l) . IStrickland & Heckmanl (120091) has also 
argued that the star-formation history of M82's starburst core 
is not well constrained before 10 Myr ago (as observed from 

3 I ignore t he relatively small difference between the up per mass limit 
of 100 M Q in Forster Schreiber et al. (2003) and 120 Mq in lCervifio eTaT] 
<200Q|) and lVoss et al.1 420091) . Since stars with masses 100 to 120 Mq can 
affect stellar diagnostics, converting to that IMF may require an adjustment 
to the star-formation history beyond a simple mass scaling. 
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FIG. 2. — History of the SLR masses in M82's "3D regi on" ISM for the star- 
formation history given in Forster Schreiber et al. (2003). The black lines are 
for 26 A l and grey lines are for 60 Fe; solid is usin g the yields in I Voss et al J 
12009) and dashes is for the Cervino et al. (2000) yields. We presently ob- 
serve M82 at t = 0; I assume there are no bursts of star-formation after then, 
so that the masses inevitably decay away. 

Earth), and ma y have extended as far back a s 60 Myr ago. 
Thus, I take the lForster Schreiber et al.l ( 120031) history merely 
as a representative example of fluctuating SFRs. 

The calculated 26 Al (black) and 60 Fe (grey) masses are plot- 
ted in Figure|2] At first, there is no SLR mass in the starburst, 
bec ause it take s a few Myr for SLR injection to start. With 
the IVoss et al l (120091) yields, the SLR masses rise quickly 
and peak ~ 5 Myr ago (as observed from Earth). The SLR 
masses drop afterwards. Yet they are still within a factor of 
1 .7 of their peak values even now, ~ 5 Myr after the last star- 
formation burst. If there is no further star-forma tion, the SLRs 
will m ostly vanish over the next 10 Myr. The ICervino et al.l 
(120001) yields predict a greater role for supernovae from lower 
mass stars, so the fluctuations are not as great; the 60 Fe mass 
remains roughly the same even 10 Myr from now. As long 
as there has been recent star-formation in the past ~ 5 Myr, 
the SLR abundances are at least half of those predicted by the 
steady-state assumption. 

There is a more fundamental reason to expect that the 
steady-state SLR abundances are roughly correct for star- 
bursts. A common way of estimating star-formati on rates in 
starbu rsts is to use the total infrared luminosity dKennicutj 
119981) . which is nearly the bolometric luminosity for these 
dust-obscured galaxies. Young stellar populations, contain- 
ing very massive stars, are brighter and contribute dispropor- 
tionately to the bolometric luminosity. Therefore, both the 
luminosity and the SLR abundances primarily trace young 
stars. To compare the bolometric luminosity, I ran a Star- 
burst99 (v6.04) model of a Z = 0.02 metallicity coeval stellar 
population with a Salpeter IMF (dN/d M oc M" 2 35 ) between 
0.1 and 120 M Q iLeitherer etaLl (119991) . I then calculate the 
SFR that would be derived from these luminosites using the 
lKennicuttl ( fl998l) conversion, and then from that, the expected 
steady-state SLR masses from equation [4] The "bolometric" 
26 Al masses are compared to the actual masses in Figure [3] 

Although the very youngest stellar populations are bright 
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FIG. 3. — How the bolometric luminosity traces 26 Al mass for a coeval stel- 
lar population with age t. The grey line is the predicted steady state 26 Al mass 
I would predict from the bolometric luminosity of the population, whereas the 
black lines are the actual mass of 26 Al (solid for[Voss et al. 2003and dashed 
for lCervino et al.120001) . 

but not yet making SLRs, the bolometric luminosity (grey) is 
a good tracer of 26 Al mass (black) for stellar populations with 
ages between 3 and 20 Myr. For most of the interval, the bolo- 
metric 26 Al masses are within a factor 2 of the act ual masses. 
For po pul ations between 15 My r and 20 Myr, the IVoss et ail 
( 120091) and lCervino e t al. ( 2000) predictions envelop the bolo- 
metric 26 Al masses. At 20 to 25 Myr old, the bolometric 26 Al 
masses are about twice the true masses. For older populations 
still, the true 26 Al masses finally die away while the bolomet- 
ric luminosity only slowly declines. Note that, if stars have 
been forming continuously for the past 100 Myr, over half of 
the luminosity comes from stars younger than 20 Myr. Thus, 
the use of the bolometric luminosities introduces a factor < 3 
systematic error. 

In short, the use of bolometric luminosity as a SFR indi- 
cator, and the natural variability in the star-formation rates of 
starbursts can lead to overestimations of the SLR abundances 
by a factor ~ 3. But I estimate the SLR abundances of true 
starbursts are a hundred times higher than in the present Milky 
Way (equation [8] and Table [TJ. The ratio is so great that the 
systematic effects do not undermine the basic conclusion that 
SLR abundances are much larger in true starbursts. 

3.2. Are SLRs mixed quickly enough into the gas? 

Although the average levels of SLRs in starbursts and high- 
z normal galaxies are high, that does not by itself mean the 
SLRs influence the environments for star-formation. While 
SLRs can play an important role in star-forming regions, by 
elevating the ionization rates and by being incorporated into 
solid bodies, SLRs trapped in ionized gas are irrelevant for 
these processes. 

The mixing of metals from young stars into the ISM gas 
mass is usually thought to be very slow in the present Milky 
Way, compared to SLR lifetimes. The massive stars respon- 
sible for making SLRs often live in star clusters, which blow 
hot and rarefied bubbles in the ISM. Supernovae also exca- 
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vate the coronal phase of the ISM dMcKee & Ostrikei|[T977h . 
Turbulence within the bubbl es mixes the SLRs and homoge- 
nizes their abundances (e.g. jMartin et al.ll2.010b . over a time 
scale f m i x w L/v mr b, where L is the outer scale of turbu- 
lence (typica l size of the largest eddies) and Vn, r h is the tur - 
bulent speed dRoy & Kuntr]ll995l : IPan & Scannapiecdl2010h . 
The large outer scale of turbulence, ^100 - 1000 pc, and 
the slow turbulent speeds (~ 5-10 km s" 1 ) in the Milky 
Way imply mixing times of ~ 10 - 200 Myr. Even if 
the 26 Al is homogenized within the superbubbles, this low 
density hot gas requires a long time to mechanically af- 
fect cold star -forming clouds, because of the large den- 
sity contrast (Ide Avillez & Mac Lowl |2002). Mixing be- 
tween the phases, particularly warm and hot gas, is acceler- 
ated by Rayleigh-Ta ylor and Kelvin-Helmholtz instabilities 
dRov & Kunthll 19951) . but o verall, mixing takes tens of M yr to 
ope rate in the Mil ky Way (Ide Avillez & Mac Lowll2002l; see 
also lClaytonll983L where mixing times are between the warm 
ISM from evaporated H I clouds, cool and large H I clouds, 
and molecular clouds). Thus, SLRs like 26 Al are thought to 
decay long before they are mixed thoroughly with the star- 
forming gas. Indeed, studies of the abundances of longer lived 
isotopes in the primordial Solar Syst em supports longer mix- 
ing times of ~ 5 - 100 Myr (e.g.. iMever & Clayton! 120001: 
Muss et al]l20"09h . 

It is thought that these obstacles existed, at least qualita- 
tively, in the z « 0.45 Milky Way, when the Solar System 
formed. These problems are part of the motivation for invok- 
ing a local sou rce of SLRs, inc luding energetic particles from 
the S un itself dLee et al.l l 19981). injection from a nearby AGB 
star (iBusso et al.1 12003]). or injection from an anomalously 
near by supernova ([Cameron & Truraiil 119771) or Wolf-Rayet 
star (lArnould et al.l (19971) . Recently, though, several authors 
proposed models that might overcome the mixing obstacle, 
where young stars are able to inject SLRs into star-forming 
clouds. A motivation behind these models is the idea that 
molecular clouds are actually int ermittent high density turbu - 
lent fluctuations in the ISM (e.g. JMac Low & Klesserill2004) . 
and the supernovae that partly drive the turbulence - indirectly 
forming the molecula r clouds - also a re the sources of SLR s 
dGounelle et al.l2009l) . In the model o flGounelle et alj d2009l) . 
old superbubbles surrounding stellar clusters form into molec- 
ular clouds after plowing through the ISM for ~10 Myr. 
Supernovae continue going off in the star clusters, adding 
their SLRs into these new born molecular clouds (see also 
IGounelle & Meynell2012l) . Simulations by IVasileiadis et afl 
d2013l) also demonstrate that SLRs from supernovae going off 
very near giant molecular clouds are mi xed thoroughly w ith 
the star-forming gas. On a different note. lPan et al.l (120121) ar- 
gued that supernovae remnants are clumpy, and that clumps 
could penetrate into molecular clouds surround star clusters 
and inject SLRs. 

If these scenarios are common, then a large fraction of the 
produced SLRs reaches the star-forming gas before decaying. 
In fact, these mechanisms may be so efficient that SLRs are 
concentrated only into star-forming molecular gas, a minority 
of the Milky Way gas mass. If so, then the abundance of SLRs 
within Galactic molecular gas (M^lr/Mui) is greater than the 
mean background level (MsLR/M gas ); in this way, SLR abun- 
dances could r each the elevated levels that existed in the ea rly 
Solar System (IGounelle et al.ll2009t IVasileiadis et al.ll2013l) . 

I note, however, that the homogeneity of the 26 Al abun- 
dance in the early Solar System is controversial; if the abun- 



dance was inhomogeneous, that is inconsistent with efficient 
SLR mixing within the Solar System's progenitor molecular 
cloud. Although IVilleneuye et all (12009 1) concl ude that 26 Al 
had a constant abundance. iMakide et alj d201 ll) find that the 
26 Al abundance varied during the earliest stages of Solar Sys- 
tem formation, when aluminum first condensed from the Solar 
nebula. 

What is even less clear, though, is how similar the Milky 
Way is to starbursts and the massive high-z normal galaxies 
that host the majority of the cosmic star formation. As in 
the Milky Way, supernovae in starbursts like M82 probably 
blast out a hot phase. But the hot phase escapes in a rapid 
hot wind in starbursts with high star-formation densities (> 
0.1 M Q yr" 1 kpc~ 2 : IChevalier & Cleggdll985b iHeckman et all 
119901) . There is direct evidence for this hot wind f rom X-ray 
observations (IStrickland & Heckmanll2007l 120091) . Further- 
more, supernova remnants are observed to expand quickly in 
M82, implying that the y are going off in a very low den- 
sity environment (e-g-> Beswick et al.l 120061; compare with 
IChevalier & Franssonll200lir Cool and warm gas is observed 
outflowing alo ngside the hotter wind, possib ly entrained by 
the hot wind dStrickland & Heckmanl 120091) . Whereas the 
edges of superbubbles eventually cool and fade back into the 
ISM in the Milky Way after a 10 - 100 Myr, any superbubble 
gas that does cool in these starbursts could still be pushed out 
by the wind. If the SLRs are trapped in the wind, they reside 
in these starbursts for only a few hundred kyr. 

But the ISM conditions in the starbursts are much differ- 
ent, with higher densities in al l phases, higher pr essures, and 
higher turbulent speeds (e.g., Smith et al. 2006). Starbursts 



are physically small, with radii of ~ 200 pc at z = - com- 
parable to the size of some individual superbubbles in the 
Milky Way. The eddy sizes therefore must be smaller and 
mixing processes could be faster. To demonstrate jus t how 
small superbubbles are in starbursts, ISilich et aD (120071) mod- 
eled the superbubble surrounding the star cluster M82 A-l in 
M82, which has a mass of 2 x 10 6 M . They find that the 
wind propagates only for a few parsecs before being shocked 
and cooled. Stellar ejecta in the core of the cluster also cool 
rapidly in their model. The turbulent mixing time is therefore 
much smaller, ~ 200 kyr for an eddy length of 10 pc and a 
turbulent speed of 50 km s" 1 . 

Conditions are even more extreme in present day compact 
ULIRG starbursts, where the ISM is so dense (~ 10 4 cm -3 ; 
iDownes & Solo mon 1998) th at a hot phase may be unable to 
form dThornton et al.ll 19981 iThompson et al.ll2005l). Instead 
the IS M is almost entirely molecular dDownes & Solomon! 
119981) . Indeed, observations of supernovae in Arp 220 in- 
dicate they are goi ng off in average density molecular gas 
dBateja t et al. 201 1). Supernovae remnants fade within a few 
tens of kyr into the ISM , due to powerful radiative losses 
dMcKee & Ostrikei|[T977h . The SLRs then are incorporated 
into the molecular ISM in a turbulent mixing time, the whole 
process taking just a few hundred kyr. The main uncer- 
tainty is then, not whether the SLRs are injected into the 
molecular gas, but whether these SLR-polluted regions of 
the molecular gas fill the entire starburst. Turbulent mixing 
smo oths abundances over regi ons the size of the largest ed- 
dies dPan & Scannapiecol2010l) . but if the distribution of SLR 
injection sites varies over larger scales, the final SLR abun- 
dance may also vary on these large scales. 

We know very little about the conditions in high red- 
shift galaxies. Star formation in main sequence galaxies 
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is dominated by massive galaxies with large star-formation 
rates at high redshift. These massive galaxies are several 
kpc in radius, but contain large amounts of molecular gas 
(iDaddi et al.ll2010h . They also have large star-formation den- 
sities and host winds. In the more extreme galaxies, radiativ e 
bremsstrahlung losses stall any hot wind (ISilich et al.ll2.010t>. 
Turbu lent speeds in these galaxies are high (IGreen et all 
120 101) . implying faster turbulent mixing than in the Milky 
Way. But it is not clear which phase the SLRs are injected 
into or how long it takes for them to mix throughout star for- 
mation regions. The effects of clustering in the ISM is also 
uncertain, but it probably is important in these galaxies, where 
huge clumps (> 10 8 M Q and a kpc wide) are observed (e.g., 
iGenzel et alJ|20lTb . 

To summarize, while there are reasons to expect that most 
of the SLRs synthesized by young stars in the Milky Way de- 
cay before reaching star-forming gas, this is not necessarily 
true in starbursts or high-z normal galaxies. Turbulent mixing 
is probably fast, at least in compact starbursts which are phys- 
ically small. On the other hand, winds might blow out SLRs 
before they reach the star-forming gas, at least in the weaker 
starbursts. Clearly, this issue deserves further study. 

4. IMPLICATIONS 

4.1. Implications for the early Solar System and Galactic 
stars of similar age 

The rapid evolution in SSFRs implies that Galactic back- 
ground SLR abundances were up to twice as high during the 
epoch of Solar System formation (4.56 Gyr ago; z ~ 0.44). If 
the evolution of the Milky Way's gas fraction was comparable 
to that in observed massive main sequence galaxies, the en- 
hancement may have been onl y ~ 50% above present values 
(see the discussion in section |2~T1 G enzel et al.l l20ldh . The 
inferred primordial abundances of 6l) Fe in the Solar System 
are in fact up to twice as high as in the cont emporary Milky 
Way, as determined with g amma-ray lines dWilliamsl 12008b 
iGounelle & Mei bom 2008). This overabundance is cited as 
evidence for an individual, rare event enriching the early So- 
lar System, or the gas that formed into the Solar System, with 
SLRs. However, my calculations show this is not necessar- 
ily the case: the twice high abundances of 60 Fe in the early 
Solar System can arise simply because the Galaxy was more 
efficient at converting gas into stars 4.5 Gyr ago. 

The primordial abundance of 26 Al in the Solar System 
was abo ut six times higher than the mean Galactic value at 
present dDiehl et al.l l2006aT) . or three times higher than the 
mean Galactic abundance at z = 0.44 assuming that equation|9] 
holds. Even so, the normal effects of galaxy evolution are a 
potential contributor to the greater 26 Al abundances, assuming 
efficient mixing of 26 Al with the molecular gas of the Milky 
Way. Furthermore, the high abundances of 26 Al in the early 
Solar System are actually typical of star-formation at z ~ 1-2 
- when most cosmic star-formation occurred. In this sense, 
the early Solar System's 26 Al abundance may be normal for 
most planetary systems in the Uni vers e. 

As I have discussed in Section 13721 it is not clear whether 
the background abundances of 26 Al and other SLRs actu- 
ally represent those of typical star-forming gas; if mixing 
takes more than a few Myr, thes e SLRs could not have 
affected star and planet formation dMeyer & Clayton] 120001 
Ide Avillez & Mac Lowll200l iHuss et alj|2009t) . But although 
there may have been a wide distribution of abundances if mix- 
ing is inefficient, the mean of the distribution is still higher 



simply because there were more supernovae and young stars 
per unit gas mass. Thus, a greater fraction of star formation 
occurred above any given threshold in the past. In addition, 
the Galactic background level can be meaningful if a large 
fraction of the SLRs fro m you n g stars make it into t he cold 
gas, and IGounelle et al.l (120091): IGounelle & Mevnetl (120121) , 
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(120121) . and lVasileiadis et alj (120131) have presented 
mechanisms where this can happen. Although there is sug- 
gestive evidence that these mechanis ms did not operate for 
the Solar System (Makide et al. 201 2), there is no reason they 
could not have worked for other star systems of similar ages. 
Then the Solar System's relatively high abundance of SLRs, 
and 60 Fe in particular, may be common for Galactic stars of 
similar ages, even if through a different process. 

Finally, as I noted, these conclusions depend on how the 
yields of SLRs from massive stars change with metallicity, 
and what the mean Galactic metallicity was at the epoch of 
Solar System formation. 

4.2. The Ionization Rate and Physical Conditions in 
Starbursts' Dense Clouds 

Radioactive decay from SLRs injects energy in the form 
of daughter particles into the ISM. The decay particles, with 
typical energies of order an MeV, ionize gas and ultimately 
heat the ISM, if they do not escape entirely. The high abun- 
dances of SLRs, including 26 Al, can alter the ionization state 
of molecular gas in these galaxies. The ionization rate, in par- 
ticular, is important in determining whether magnetic fields 
can thread through the densest gas. I focus here on the con- 
tribution from 26 Al, which dominated the ionization rate from 
SLRs in the primordial Solar System (Um ebavashi & Nakanol 
120091) . For the sake of discussion, I assume that the SLRs 
are well -mixed into the gas, despite the uncertainties (sec- 
tion E2). 

Each 26 Al decay releases an energy £decay i nto tne ISM 
in the form of MeV positrons and gamma rays. If each 
atom in the ISM takes an energy E^ on to be ionized, each 
26 Al decay can therefore ionize -Edecay/^ion atoms. In 82% 
of the decays, a positron with kinetic energy of 1.16 MeV 
is released, and the po sitron is slowed by ionization losses 
dFinocchi & Gail|[T997h . The minimum energy per decay, af- 
ter accounting for this branching ratio, that goes into ioniza- 
tion is Zs™" = 0.95 MeV, when the medium stops the positron 
(inflight annihilation los ses are negligible at these energies; 
iBeacom & Yiikse l 2006). The annihilation of the positron 
with electron in the ISM produces gamma rays of total en- 
ergy 2 x 0.511 = 1.022 MeV. In addition, in very nearly all 
26 Al decays, a 1.809 MeV gamma ray is produced. These 
gamma rays only interact with the ISM over columns of sev- 
eral g cm" 2 , so only in parti cularly dense regions w ill they 
contribute to the ionization dFinocchi & Gaill|1997l) . When 
they do, E^ ay = 3.60 MeV. iStepinskj d 19921) gives E ioa as 
36.3 eV, so that the ionization rate is 



c 



•^Al-26^d 



Al-26 = 



ecay 



(36.3 eV)T deC ay 



(12) 



In terms of gas consumption time, the ionization rate from 
26 Al is 



C = (1.4-5.1) x 10" 18 sec" 



20 Myr 



(13) 



My results for the mean ionization rate from 26 Al of some 
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characteristic starbursts are shown in Table Q~] they are in 
the range 10~ 18 - 10~ 17 sec" 1 . The maximal ionization rates 
are roughly an order of magnitude higher than that found in 
early Solar System, a dense environment with £ ps (0.6- 1) x 
IP" 18 sec -1 dFinocchi & Gaillll997t lUmebavashi & Nakanol 
l2009h . 

Even if the 26 Al is in the cold star-forming gas, it could 
actually be condensed into dust grains instead of existing in 
the gas phase. Yet the decay products still escape into the ISM 
from within the grain. The attenuation of gamma rays at ~ 1 
MeV is dominated by Compton scattering, requiring columns 
of a few g cm" 2 to be important. Thus, gamma rays pass freely 
through dust grains that are smaller than a centimeter. The 
energy loss rate of relativistic electrons or positrons in neutral 
matter is approximately 

dK 9 , v -> p;Zj I" K+m e c 2 2, m e c 2 ~\ ,„ JN 
— = -m e c 2 a T > ■^ L - L In ^— +-ln— — (14) 

ds 4 Ajm H [ m e c 2 3 (Ej) \ 

from lSchlickeiserl (120021) . where K is the particle kinetic en- 
ergy, s is the physical length, and aj is the Thomson cross 
section. The sum is over elements j; for heavy elements Z, ps 
Aj/2, pj is the partial density of each element within the grain, 
and (Ej) is related to the atomic properties of the element. I 
take the bracketed term to have a value ~ 5 and pw3g cm" 2 . 
Then the stopping length is K/(dK/ds) ps 0.3 cm (K/MeV), 
much bigger than the typical grain radius. Thus, 26 Al (and 
other SLRs) in dust grains still contribute to the ionization of 
the ISM. 

On the other hand, are the positrons actually stopped in 
starburst molecular clouds, or do they escape? For neu- 
tral interstellar gas, the stopping column of MeV positrons 
is Kp/(dK/ds) ps 0.2 g cm" 2 ( ISchlickeiserj|2002l) . The col- 
umn densities of starb ursts range from ~ 0.1 - 10 g cm" 2 
(e.g., iKennicutl 119981) . and the columns through individ- 
ual molecula r clouds are expecte d to be similar to those of 
the galaxies ([Hopkins et al.ll20l2h . In the denser molecular 
clouds, positrons are stopped even if they are not confined 
at all. In ma ssive main sequen ce galaxies, the columns are 
~ 0. 1 g cm" 2 dDaddi et al.l2010l) . insufficient to stop positrons 
moving in straight lines. If magnetic waves with a wavelength 
near the positron gyroradius scale exist in these clouds, they 
efficiently scatter positrons and confine them. However, these 
waves are probably d amped quickly in dense neutral gas (see 
iPrantzos et al.l l201ll and references therein). On the other 
hand, positrons move along magnetic field lines, and if the 
lines themselves are twisted on a scale Ag, the positrons are 
forced to random walk with a similar mean free path. As long 
as Xb is less than about a third of the molecular cloud size, 
then positrons are stopped in these galaxies. 

If it is well mixed with the molecular gas, does 26 Al dom- 
inate the ionization rate in molecular gas in starbursts, and 
if so, what physical conditions does it induce? The star- 
burst 26 Al ionization rates are about an order of magnitude 
lower than the canonical cosmic ray-sustained ionization rates 
in most Milky Way molecular clouds, but in some of the 
densest Galactic starless cores, the ionization rate drops to 
- 10" 18 sec" 1 (ICaselli et al.l l200l iBerein & Tafallal [2007b . 
Cosmic rays in starbursts can sustain much higher ionization 
rates (up to ~ 10~ 14 sec" 1 ; c.f.. lPapadopoulosl2010l) . but cos- 
mic rays can be deflected by magnetic fields, possibly pre - 
venting them from penetrating high columns. In lLackil (120 1 2|) . 



I found that the gamma-ray ionization rate in starbursts can 
therefore be anywhere from 10~ 22 - 10~ 16 sec" 1 , with values 
of- (1-3) x 10" 19 sec" 1 inM82and-(5-8) x 10" 17 sec" 1 
in Arp 220's radio nuclei. In the dense clouds of most star- 
bursts, 26 Al radioactivity could exceed the ionization rate 
over gamma rays, setting a floor to the ionization rate. In 
the most extreme starbursts, with me an gas surfac e densities 
of > 3 g cm" 2 (c.f. equation 10 of lLackil 120 121) . however, 
gamma-ray ionization is more important, since the gamma- 
ray ionization rate depends strongly on the density of gas and 
compactness of the starbursts. Unlike the uncertainty of how 
SLRs and their positron decay products are transported and 
mixed with the gas of starbursts, gamma rays propagate rel- 
atively simply, so the gamma ray ionization rates are more 
secure. 

An 26 Al-dominated ionization rate has implications for the 
physica l cond itions of star-forming clouds. According to 
iMcKed ( 119891) . the ionization fraction of a cloud with hydro- 
gen number density n# is 

* e ps 1.4 x 10" 8 f tS r) ( z ) (15) 

V 10" 18 sec ~ l J V10 4 cm" 3 / v ; 

when the ionization rate is low. We see that the ionization 
fraction of cloud with density kh = 10 4 cm" 3 are (1-4) x 10~ 8 , 
if the ionization is powered solely by 26 Al decay. For these 
ionization fractions, the ambipolar diffusion time of a molec- 
ular core, the time for magnetic fields to slip from the gas, is a 
few times its free-fall time. Since clouds with strong magnetic 
fields do not collapse until the field slips away by ambipo- 
lar di ffusion dMestel & SpitzedfT956T: iMouschovias & Spitzed 
fl976l) . this means that 26 Al-ionized clouds in starbursts col- 
lapse quasi-statically, as in the Milky Way. 

On the other hand, the energy injectio n from 26 Al has essen- 
tially no effect on the gas temperature. iPapadopoulosI (120101) 
gives the minimum gas temperature of gas as: 

7T n = 6.3 K [(0.0707n4 /2 C-i8 + 0.186 2 n^) 1/2 -0.186«4 /2 ] 2/3 , 

(16) 

which was derived under the assumption that there is no heat- 
ing from interactions with dust grains or the dissipation of 
turbulence in the gas, for gas with density n\ = (n# /10 4 cm" 3 ) 
and ionization rate £_ig = (£/10~ 18 sec" 1 ). In typical star- 
bursts, I find that 26 Al decay alone heats gas of density = 
10 4 cm " 3 to 7 2-5 K (0.1-0.5 K for n H = 10 6 cm" 3 ). As I 
note in lLackil (120121) . under such conditions, dust heating is 
more likely to set the temperature of the gas than ionization, 
raising the temperature to the dust temperature for densities 
> 40000 cm" 3 . 

5. CONCLUSIONS 

The high SSFRs of starbursts and high-z normal galax- 
ies implies high abundances of 26 Al and other SLRs in their 
ISMs. In true starbursts, these abundances are enormous, with 
X( 26 A1) ps 10" 9 and 26 A1/ 27 A1 ps 10" 3 . The SSFRs of normal 
galaxies evolve rapidly with z; even taking into account higher 
gas fractions, the SLR abundances were about 3-10 times 
higher at z ~ 2 than in the present Milky Way. Even at the 
epoch of Solar System formation, the mean SLR abundance 
of the Milky Way was 1 .5 to 2 times as high as at the present 
(Figure [TJ. This alone could explain the high abundances of 
60 Fe in the early Solar System, and reduce the discrepancy in 
the 26 Al abundances from a factor of ~ 6 to — 3. In this way, 
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the cosmic evolution of star-forming galaxies may have di- 
rect implications for the early geological history of the Solar 
System. The first main uncertainty is whether the SLRs pro- 
duced by massive stars is well-mixed with the molecular gas: 
they may instead be ejected by the galactic winds known to be 
present in starbursts and high-z galaxies, or decay before they 
can propagate far from its injection site, or before it penet rates 
cold gas. I discussed these uncertainties in section 13.21 The 
other uncertainty is how SLR yields depend on metallicity. 

The most direct way to test the high 26 Al abundances of 
starburst galaxies is to detect the 1.809 MeV gamma-ray line 
produced by the decay of 26 Al, directly informing us of the 
equilibrium mass of 26 Al. Whether most of the 26 Al is ejected 
by the superwind can be resolved with spectral informa- 
tion. The turbulent velocities of molecular gas in starbursts, 
~ 100 km s" 1 dDownes & Solomon! 1 19981) . is much smaller 
than the bulk speed of the superwind, which is hundreds or 
even thousands of km s" 1 ( IChevalier & Clegglll985l) . Unfor- 
tunately, the 26 Al line fluxes predicted for even the nearest 
external starbursts (~ 10~ 8 cm " 2 sec" 1 ) are too l ow to detect 
with any planned in strument ( lLacki et al.l 120121) . However, 
ICrocker et al.l (1201 ll) have argued that the inner 100 pc of the 
Galactic Center are an analog of starburst galaxies, launch- 
ing a strong superwind. The 26 Al line from this region should 
have a flux of ~ 2 x 10~ 5 cm" 2 sec" 1 , easily achievable with 
possible futur e MeV telesc opes like Advanced Compton Tele- 
scope (ACT; iBoggsl l2006h or Gamma-Ray B urst Investiga- 
tion v ia Polarimetry and Spectroscopy (GRIPS : lGreiner et al.l 
1201 l!) . A search for the 26 Al signal from this region may in- 
form us on its propagation, since it is nearby and resolved 
spatially If the 26 Al generated in the Center region is actu- 
ally advected away by the wind, the "missing" 26 Al should 
be visible several hundred pc above and below the Galactic 
Plane near the Galactic Center Q. Resolved measurements of 
the Galactic Center can also inform us on whether the 26 Al is 
present in all molecular clouds in the region (and therefore is 
well-mixed), or just if it is trapped near a few injection sites. 

If the SLRs do mix with the star-forming molecular gas 
of these galaxies, there are implications for both their star 
formation and planet formation. Ionization rates of 10~ 18 - 



10 17 sec l , like those in some Milky Way starless cores, 
result from the energy injection of 26 Al decay in starbursts. 
While cosmic ray ionization rates can easily exceed those ion- 
ization rates by orders of magnitude in gas they penetrate into, 
and while gamma rays produce higher ionization rates in the 
most extreme starbursts like Arp 220, 26 Al might dominate the 
ionization rate in the dense clouds shielded from cosmic rays 
in typical starbursts. In starbursts' protoplanetary disks, 26 Al 
can provide moderate ionization through all columns, possi- 
bly elimi nating the "dead zones" where ther e is little accre- 
tion (e.g.. lGammiel fT996; Fatuzzo et al. 2006). Any planetes- 
imals that do form in starbursts may have much higher radio- 
genic heating from 26 Al. Admittedly, studying the geological 
history of planets, much less planetesimals, in other galaxies 
is very difficult for the forseeable future. However, at z ~ 2 
the Milky Way likely had background SLR abundances ^10 
times higher than at present, so the effects of elevated SLR 
abundances may be studied in planetary systems around old 
Galactic stars. 

On a final point. iGilmour & Middletonl (120091) propose that 
the elevated abundances of 26 Al in the early Solar System are 
mandated by anthropic selection, since high SLR abundances 
are necessary for planetesimal differentiation and the loss of 
volatiles, but that explanation may be difficult to maintain. 
If high 26 Al abundances, far from being very rare, are actu- 
ally typical of high-z and starburst solar systems (and indeed, 
much of the star-formation in the Universe's history), the an- 
thropic principle would imply that most technological civi- 
lizations would develop in solar systems formed in these en- 
vironments. Instead of asking why we find ourselves in a sys- 
tem with an 26 Al abundance just right to power differentiation 
and evaporate volatiles, we must ask why we find ourselves in 
one of the rare solar systems with sufficient 26 Al that formed 
in a normal spiral galaxy at z ~ 0.4, instead of the common 
26 Al-enriched solar systems formed at z ~ 2 or in starbursts. 
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